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O ' Abstract 

cN : 

, We prove a fractional Noether's theorem for fractional Lagrangian systems invariant 

j^ ' under a symmetry group both in the continuous and discrete cases. This provides an 

explicit conservation law (first integral) given by a closed formula which can be algorith- 

mically implemented. In the discrete case, the conservation law is moreover computable 

>0 \ in a finite number of steps. 
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c^ ■ 1 Introduction 

A conservation law for a continuous or a discrete dynamical system is a function which is 
constant on each solution. Precisely, let us denote by g = (5(i))te[a,fe] (resp. Q = iQk)k=o,...,N) 
^ , the solutions of a continuous (resp. discrete) system taking values in a set denoted by E. 

^^ I Then, a function C : i? — t- M is a conservation law (or first integral) if, for any solution q 

pi ' (resp. Q), it exists a real constant c such that: 

^ C{q{t))=c,yte[a,b]. (resp. C(Qfc) = c, VA; = 0,...,iV.) (1.1) 

o' 

CS| ' Conservation laws are generally associated to some physical quantities like total energy or 

angular momentum in mechanical systems and sometimes, they can be used to reduce or 
integrate by quadrature the equation. 

X 

^ . In this paper, we study the existence of classical conservation laws for continuous or discrete 

fractional Lagrangian systems introduced in [1] and [1] respectively. We follow the usual 
strategy of the Noether's theorem which provides explicit conservation laws for classical La- 
grangian systems invariant under symmetries. In particular, we prove a fractional Noether's 
theorem both in the continuous and discrete cases providing explicit formula of the conser- 
vation law. This formula is algorithmic so that arbitrary hight order approximations can be 
computed. In the discrete case, the algorithm is moreover finite. 

Previous results in this direction have been obtained by several authors ([3], [5], [6]). However, 
in each of these papers, the conservation law is not explicitly derived in term of the symmetry 
group and the Lagrangian but implicitly defined by functional relations ([5], [6]) or integral 
relations [3]. 



The paper is organized as follows. In Section [2l we first remind classical definitions and 
results concerning fractional derivatives and fractional Lagrangian systems. Then, we prove 
a transfer formula leading to the statement of the fractional Noether's theorem based on a 
result of Cresson in [2] . Finally, we give some remarks concerning the previous studies on the 
subject and we give some precisions on our result. In Section [3l we first remind the definition 
of the discrete fractional Lagrangian systems in the sense of [1] . Then, introducing the notion 
of discrete symmetry for these systems, we state the discrete fractional Noether's theorem. 
We conclude with some remarks and a numerical implementation. 

2 A fractional Noether's theorem 

In this paper, we consider fractional differential systems in M where (i S N* is the dimension. 
The trajectories of these systems are curves q in C^{[a^ h\^W^) where a < b are two reals. Let 
< a < 1 denote a fractional order. 

2.1 Reminder about fractional Lagrangian systems 

We first review classical definitions extracted from (8l \T0\ \TT\ of the fractional derivatives 
of Riemann-Liouville. Let T be the Gamma function of Euler and / be an element of 
C^{[a,b],]R. ). The fractional left (respectively right) integral of order a with inferior limit 
a (respectively superior limit b) of / is given by: 

VtG]a,6], r/(t) = -l^ f\t-yr-'f{y)dy, (2.1) 

r(a) Ja 

respectively: 

Vt G [a, b[, lXf{t) = -i^ j\y - tr-'f{y)dy. (2.2) 

r(a) Jt 

Then, the fractional left (respectively right) derivative of order a with inferior limit a (re- 
spectively superior limit b) is given by: 

ViG]G,6], D^/(i) = I (/!-"/) (i), (2.3) 

respectively: 

Vt € [a,bl Dlf{t) = -A(^4-"/)(t). (2.4) 

A fractional Lagrangian functional is an application defined by: 

£" : C2([a,6],M°') -^ M (2.5) 

q ^ / L{q{t),Dlq{t),t)dt, 

J a 

where L is a Lagrangian i. e. a C^ application defined by: 

L: W^ xW^ y.[a,b] — > M (2.6) 

{x,v,t) I — > L{x,v,t). 

It is well-known that extremals of a fractional Lagrangian functional can be characterized by 
a variational principle as solutions of a fractional differential system, see [1]. Precisely, g is a 



critical point of £° if and only if q is solution of the fractional Euler-Lagrange equation given 
by: 



r)T f BT \ 



A dynamical system governed by an equation shaped as ()EL"P is called fractional Lagrangian 
system. 

For a = 1, we have -Di = — -D+ = d/dt and ()EL"P coincides with the classical Euler-Lagrange 
equation (see [2]). 

2.2 Symmetries 

We first review the definition of a one parameter group of diffeomorphisms: 

Definition 1. For any real s, let <p{s, •) : M — > M be a diffeomorphism. Then, $ = 
{<j){s, •)}sgR is a one parameter group of diffeomorphisms if it satisfies 

1. 0(0,-) = Idf,d, 

2. Vs, s' € M, 0(s, •) o </.(s', •) = ^{s + s', •) , 

3. (/) is of class C^ . 

Classical examples of one parameter groups of diffeomorphisms are given by translations and 
rotations. 

The action of a one parameter group of diffeomorphisms on a Lagrangian allows to define the 
notion of a symmetry for a fractional Lagrangian system: 

Definition 2. Let ^ = {(j){s, •)}s£R be a one parameter group of diffeomorphisms and let L 
he a Lagrangian. L is said to he W^-invariant under the action of ^ if it satisfies: 

^q solution of (IEL^ . Vs G M, LU{s,q),D'l{cp{s,q)),t\ =L{q,D^q,t). (2.7) 

2.3 A fractional Noether's theorem 

Cresson [5] and Torres et al. [6] proved the following result: 

Lemma 3. Let L he a Lagrangian Vi^ -invariant under the action of a one parameter group 
of diffeomorphisms <l> = {<j){s, ■)}s£R- Then, the following equality holds for any solutions q 
of dEE]); 



In the classical case a = 1 , the classical Leibniz formula allows to rewrite (j2.8p as the derivative 
of a product. Precisely, from Lemma [3l it leads to the classical Noether's theorem given by: 



Theorem 4 (Classical Noether's theorem). Let L be a Lagrangian d/dt-invariant under the 
action of a one parameter group of diffeomorphisms ^ = {(/>(s, •)}sgiR. Then, the following 
equality holds for any solutions q of (EL^); 



d ( dd) , , dL , , , 



(2.9) 



where q is the classical derivative of q, i.e. dq/dt. 



Theorem S] provides an explicit constant of motion for any classical Lagrangian systems ad- 
mitting a symmetry. In the fractional case, such a simple formula allowing to rewrite (|2.8|) as 
a total derivative with respect to t is not known. The next theorem provides such a formula 
in an explicit form: 



Theorem 5 (Transfer formula). Let f,g €z C°°([a,6],M ) assuming the following condition 
(C): 

the sequences of functions {I^~"f ■ g^^^)p£f^* and 
(j(p) ■ I^"g)p^^, converge uniformly to on [a,b]. 

Then, the following equality holds: 



Dt!-,-!D%, = - 



y^(_iy7'^'+l-" f . qW _)_ fW . jr+l-a^ 



.r=0 



(2.10) 



Proof. Let f,gG C°^{[a, b],M. ). By induction and using the classical Leibniz formula at each 
step, we prove that for any p G N*, the following equalities both hold: 



D^f.g={-irF-y.9^P^ + 



d_ 
di 



■p-i 



Y^i-iyr+^-f.gir) 



r=0 



and 



f-Dlg = /(^) • /r"<7 + I 



■p-i 

Y^fir).jr^+l~ag 
r=0 



(2.11) 



(2.12) 



For any p € N*, we denote by Up the following function: 



p-i 



jjp .- 



^{-lyr+'-'^f-g 



• nM. 



r=0 



(2.13) 



According to ()2.1ip . for any p € N*, iip = D"_f • g — (—1)^/^ "/ • g'^ • Hence, the assumption 
made on the sequence of functions (^lP^°'f ■ g(P^)p^fq* implies that the sequence ('Up)p(=N* 
converges uniformly to DZf ■ g on [a,b]. Moreover, let us note that (lip)pGN* point- wise 
converges in t = a. Finally, we can conclude that the sequence (lip)pgN* converges uniformly 
on [a, b] to a function u equal to 



■.= ^(-iYr+^--f.g 



.nM 



(2.14) 



r=0 



and satisfying u = D'^f ■ g. Similarly, one can prove that: 



di 



J2f{r).jr^+l~ag 



.r=0 



-f-D^g. 



D 



Thus, combining Lemma [3] and this transfer formula, we prove: 

Theorem 6 (A fractional Noether's theorem). Let L be a Lagrangian D'^-invariant under 
the action of a one parameter group of diffeomorphisms $ = {<j){s, •)}seR- Let q be a solution 
of ()EL"P and let f and g denote: 



f=^{0,q) and g = ^{q,D'^q,t). 
If f and g satisfy the condition (C), then the following equality holds: 



d_ 
di 



y^ t-iy J'T+^-o^ f . g(r) _|_ f(r) . J^+l-"^ 



,r=0 



0. 



(2.15) 



(2.16) 



This theorem provides an explicit algorithmic way to compute a constant of motion for any 
fractional Lagrangian systems admitting a symmetry. An arbitrary closed approximation of 
this quantity can be obtained with a truncature of the infinite sum. 

Previous results in this direction have been obtained by several authors. We give some details 
and comments in the following. 

• In [6], Torres et al. define a fractional conservation law for fractional Lagrangian sys- 
tems. Precisely, denoting by Dq the bilinear operator 



V/,5 e C\[a,b],R''), D^{f,g) = D'^f ■ g - f ■ D^g, 
they give the following definition: a function C : W^ x 



(2.17) 



^'^^ X M*^ X [a, b] — > R 

{x,v,t) I — > C{x,v,t) 
is a fractional- conserved quantity of a fractional system if it is possible to write C in 
the form C{x,v,t) = C^{x,v,t) ■ C^{x,v,t) with D^{C^{q,D'^q,t),C^{q,D'^q,t)) = 
for any solution q of the considered system. They deduce easily that d(j)/ds{0,x) ■ 
dL/dv{x,v,t) is a fractional-conserved quantity of ()EL"p . Two important difficulties 
appear: 

— Although all these notions and results coincide with the classical ones when a = 1, 
we do not have that Dq (/,(/) = implies f ■ g is constant. 

— The decomposition C = C^ • C^ is not unique and it leads to many issues. Indeed, 
the scalar product is commutative and then C = C^ ■ C^ = C^ ■ C^ . Nevertheless, 
the operator Dq is not symmetric! Hence, in the result of Torres et al, we have: 



^, . dcj) . . dL dL ^ d(j) . . 

C{x,v,t) = ^(0,x) • -g^{x,v,t) = -g^{x,v,t) ■ ■^(0,x) 

is a fractional-conserved quantity of ()EL"P but we do not have: 



^o{f(?,^-^,t),|^(0,, 



0. 



(2.18) 



(2.19) 



• In [3], Atanackovic et al. were interested by symmetries modifying also the time variable 
t. They finally obtained in [31 Theorem 15] a formulation of constants of motion for 
fractional Lagrangian systems admitting such symmetries. Nevertheless, if we consider 
a symmetry $ which does not modify the time variable, the result of Atanackovic 
et al. is the following: for any solutions q of ()EL°p . let / = d(f)/ds{0,q) and g = 
dL/dv{q,D'i^q,t), the following element t i — > f D'^f ■ g — f ■ D'^gdy is constant on 
[a, 6]. This result is then unsatisfactory because the constant of motion is not explicit 
in contrary to our result. 



2.4 Concerning the condition (C) 

Condition (C) could seem strong or too particular. In order to make it more concrete and 
understand what classes of functions satisfy (C), we give the two following sufficient condi- 
tions: 

Proposition 7. Let f,ge C°°{[a,b],M.'^). 

1. If f and g satisfy the two following conditions: 
{h-t)P"^ 



max 

iG[a,b] 



-l/^^HOI 



and 



max 

te[a,b] 



p-1 



{t - a) 
{p-l)\ 



5(^)(t)l 



p— >-oo 



(2.20) 



{p — 1)! ■" ■ y p^oo 
then f and g satisfy the condition (C). 
2. If there exists Af > such that: 

Vp G N*, Vt € [a, 6], \f^'P\t)\ < M and \g^P\t)\ < M 
then f and g satisfy the condition (C). 
Proof. 1. Let us prove that the sequence of functions (^I^'^f ■ g^'P')p^fqt converges uniformly 



(2.21) 



to 0. The proof is similar for the sequence (/'^^ -/^ "(7)pgN*- Let us denote M = max(|/(t)|). 

t£[a,b] 

For any p € N* and any t € [a, b]: 
I'rV{t)g^\t) 



/W(t) 



T{p 



Q 



(t-y) 



p~a~l 



f{y)dy 



< M: 



(t-a) 



p—a 



(2.22) 



r(p+l 



a) 



g(p)it)\. 



Finally, as T{p + l-a)>{p- 1)! r(2 - a): 



irV{t)g(p){t] 



which concludes the proof. 

2. It follows from the first result. 



< M 



< M 



{t - a)i-" (t - a)P-i 



r(2 



a) 



{p-l)\_ 



(6-a)i-"(t-a)P-i 



r(2 



a] 



[p-l)\ 



g(p)(t)| 



(2.23) 



D 



For example, if / is polynomial and g is the exponential function, one can prove that / and 
g satisfy the point 2 of Proposition [7] and then the condition (C). If g{t) = 1/t on [a, 6] with 
a > 0, one can notice that the point 2 of Proposition [7] is not satisfied anymore. Nevertheless, 
the point 1 is true and consequently the functions / and g satisfy the condition (C). In general, 
every couple of analytic functions satisfy condition (C) and consequently ()2.10p is true. 



3 A discrete fractional Noether's theorem 

We study the existence of discrete conservation laws for discrete fractional Lagrangian systems 
in the sense of [1]. Using the same strategy as in the continuous case, we introduce the notion 
of discrete symmetry and prove a discrete fractional Noether's theorem providing an explicit 
computable discrete constant of motion. 

3.1 Reminder about discrete fractional Lagrangian systems 

We follow the definition of discrete fractional Lagrangian systems given in [3] to which we 
refer for more details. 

We denote hy N £ N*, hy h = {b — a)/N the step size of the discretization and by r = 
{tk)k=Q,...,N the following partition of [a,b]: 

\/k = 0, ...,N, tk = a + kh. (3.1) 

Let us denote by A" and A" the following discrete analogous of W^ and -D" respectively: 

(3.2) 



A^ 


: {W^Y+^ 




Q 


A^: 


(^^d-^N+l 




Q 



r,d\N 



V ''=0 / fc=l...,Af 



N 
N-k 



(3.3) 



■7- Yl O^rQk+r 



'-=0 / k=O....N-l 



where the elements (ar)reN are defined by ao = 1 and 

Vr e N*, ar = (-o^)i^ --)■■■(- -^--) _ (3.4) 

These discrete fractional operators are approximations of the continuous ones. We refer to 
|10j for more details. 

A discrete fractional Euler-Lagrange equation, of unknown Q G (M ) ~''^, is defined by: 

— (g,A-Q,r) + A:^(^— (Q,A^Q,r)j =0, (EL^) 

where L is a Lagrangian. In such a case, we speak of a discrete Lagrangian system. 



The terminology is justified by the fact that solutions of (EL? ) correspond to discrete critical 



points of the discrete fractional Lagrangian functional defined by: 

CI : (M'^)^+i -^ M (3.5) 

N 

Q ^ /iJ^L(Qfc,(A-Q)fc,tfc). 



k=l 



We refer to [1] for a proof. 

In the case a = 1, the discrete operators Aj_ correspond to the imphcit and exphcit Euler 
approximations of d/dt and (EL^) is just the discrete Euler-Lagrange equation obtained in 
019]. 

3.2 Discrete symmetries 

Here again, a discrete symmetry of a discrete fractional Lagrangian system is based on the 
action of a one parameter group of transformations on the associated Lagrangian: 

Definition 8. Let $ = {(/>(s, •)}sgR he a one parameter group of dijjeomorphisms and let L 
be a Lagrangian. L is said to be A'i^-invariant under the action of ^ if it satisfies: 



\/Q solution of pL^ , Vs G M, LU(s, Q), A° ((/)(s, Q)), rj =L(Q,A°Q,t). (3.6) 

Then, we can prove the following discrete version of Lemma [3l 

Lemma 9. Let L be a Lagrangian A'^-invariant under the action of a one parameter group 
of diffeomorphisms ^ = {(p{s, OIsgm- Then, the following equality holds for any solutions Q 

of ( ra .- 

A° (1^(0, g)) • ^{Q, A^Q, r) - |^(0, Q) • A^ (^^{Q, A^Q, r)) = 0. (3.7) 

Proof. This proof is a direct adaptation to the discrete case of the proof of Lemma [31 Let 
Q € (M°')^"'~^ be a solution of (|EL?h. Let us differentiate equation ()3.6p with respect to s and 



invert the operators A° and d/ds. We finally obtain for any s G M and any k € {1,...,A^ — 1}: 
dL 



^-'S<'-«/, av 



(t)is,Qk),A'^{(t>{s,Qk)),tk 



dL 



+ TT ('^(^' ^'^)' ^- ('^(^' ^'^)) ' ^i') ■ 7r(«' ^'^) = 0- (3 



5rE V 



9s 



Since 0(0, •) = /^Kd and Q is solution of (EL^ ), taking s = in ()3.8p leads to ()3.7p . 
3.3 A discrete fractional Noether's theorem 



D 



Let us remind that our aim is to express explicitly a discrete constant of motion for discrete 
fractional Lagrangian systems admitting a discrete symmetry. As in the continuous case, our 
result is based on Lemma [9j Let us note that the following implication holds: 



VF G M^+\ AiF = ^ 3c G M, V/c = 0, ..., N, F^ = c. 



(3.9) 



Namely, if the discrete derivative of F vanishes, then F is constant. Our aim is then to write 
()3.7p as a discrete derivative {i.e. as AL of an explicit quantity). 



We first introduce some notations and definitions. The shift operator denoted by a is defined 
by 

a: (M'^)^+i -^ (]R'^)^+i (3.10) 

Q I > cr{Q) = {Qk+l)k=0,...,N 



with the convention Qn+i = 0. We also introduce the following square matrices of length 
{N + 1). First, Ai = —Idj\[+i and then, for any r G {2,...,A^ — 1}, the square matrices 
Ar G J^N+i defined by: 



r 



yi,j = 0,...,N, (Ar) 



«J 



. {Ar)N-l,j 

where 6 is the Kronecker symbol. 

For example, for N = 5, the matrices A^ are given by: Ai 



^{j=0}^{r<i} - ^{0<i-j<r-l}^{l<j<N-r} if I <i < N -1 



if i = 
if 1 <i 

if i = N 



(3.11) 



-I dp. and 



/ 











\ 





1 ° 










\ 






-1 










-1 





1 
1 


-1 




-1 
-1 




-1 






,A3 = 




1 


-1 
-1 


-1 
-1 






1 








-1 







1 





-1 





V 1 








-1 


0/ 




V 1 





-1 


0/ 



,A4: 



/O 0000\ 

'0-10000* 

0-10000 

0-10000 

1-10 

\1 -10000/ 



Considering these previous elements, we state the following result: 

Theorem 10 (A discrete fractional Noether's theorem). Let L he a Lagrangian A":-invariant 
under the action of a one parameter group of diffeomorphisms $ = {4>{s, •)}sgiK. Then, the 



following equality holds for any solutions Q of (EL^).- 

'N-l 



Ai 



j:arAJ^iO,Q)-a\^iQ,A-_Q,r))') 

r=l ^ ^ 



0. 



(3.12) 



Combining (j3.9p and ()3.12|) , Theorem [10] provides a constant of motion for discrete fractional 
Lagrangian systems admitting a symmetry. Let us note that the discrete conservation law is 
not only explicit but computable in finite time. For an example, we give a numerical test in 
the next section. 



Before giving the proof, we remark that Theorem [TU] takes a particular simple expression 
when a = 1. Indeed, since ar = for any r > 2, ai = —1 and Ai = —IdN+i in this case, we 
obtain: 

Theorem 11 (Discrete classical Noether's theorem). Let L he a Lagrangian A\-invariant 
under the action of a one parameter group of diffeomorphisms <1> = {0(s, •)}sgiK. Then, the 
following equality holds for any solutions Q of (EL^); 



Al_ 



f(o,g)-a(f(g,AiQ,r: 



0. 



(3.13) 



This result is a reformulation of the discrete Noether's theorem proved in [3 [9]. It corresponds 
to our Lemma [9] with q = 1 using the following discrete Leibniz formula: 



VF, G G (M"')^+\ Al {F • a{G)) = AIF -G-F- A\G. 



(3.14) 



Proof of Theorem 1 1 01 According to Lemma [H equation ()3.7p holds for any A; = 1,...,A^ — 1. 
Let F = d(j)/ds{0,Q) and G = dL/dv{Q,A'^Q,T). Let us multiply ([HZD by /i" and obtain 
the following equality for any A; = 1,...,A^ — 1: 



^N-k 



\r=0 



. r=0 



[ Y, O^rFk-r ]-Gk-Fk-[Y^ ^rGk+r ] = 0, 

and as ao = 1: 

ai{Fk-l -Gk- Fk- Gk+l)k + ^ OirFk-r ' Gfc - -Ffc • ^ O-rGk+r 1 = 0, 



(3.15) 



/iV-fc 



(3.16) 



vr=2 



then: 



1 



a,A\{F-a{G))^ = - 



h 



. r=2 



^N-k 



"YarFk-rj -Gk- Fk- i^^arGk+rj ■ (3.17) 

Vr=2 / \r=2 /. 

Let us denote Jk the right term of (j3.17p for any A; = 1,...,A^ — 1. We are going to write Jk 
as the discrete derivative of its discrete anti-derivative: it corresponds to the discrete version 
of the method of Atanackovic in [3]. For any k = 1, ..., A^ — 1, we obtain Jk = {A^H)k where 
Hi := /iX^j^i Jj for any i = 0, ...,N with the convention Hq = J/v = 0. Let us provide an 
explicit formulation of the element H = (i:^j)j=o,...,Ar. 
• Case l<i < N -1: 






N^j 



y CX^ J7 'j — If ' \Sf -I / (JCrf JT a ' \_Z •I'-L-i' 



r=2 



■t ] 



N-j 



As Ej=2 Er=2 = Er=2 Ej=r' ^c have: 



j=2 r=2 j=l r=2 

(3.18) 



i 3 



Y,Y.''rFj~r-G,=Y.Y.''rF3-r-G,=Y.Y.^rFyG,^r = Y.Y.^r{F.a^{G))^. 
j=2 r=2 r=2 j=r r=2 j=0 r=2 j=0 

(3.19) 

±iien, db Z^j=i Z^r=2 ~ ^i=l Z^r=2 "■" Z^j=l Z^r=N+l-i ~ l^r=2 Z-^j=l "•" Z^r=N+l-i Z^j=l ' 
i i-r N-i i N-l N~r 

Hi = Y.Y.''r{F ■ a^{G))^ -Y,Y."r{F ■ a^G))^ - Y. Y ""ri^ ■ ct^{G)) ^. (3.20) 

r=2 j=0 r=2 j=l r=N+l-i j=l 



Finally, we have: 



Af-l N 

Hr=YY.'^rAr{i,j){F.a''{G))^, 

r=2 j=0 



(3.21) 



where the elements (^Ar{i,j)) are defined for r = 2,...,N — 1 and j = 0,...,N as the real 
coefficients in front of ar [F ■ a^ (G)) ,• • Our aim is then to express the values of these elements. 
From (|3.20p . we have for any r = 2, ..., A — 1 and any j = 0, ..., A^: 



Ar{i,j) — S{r<i}^{0<j<i-r} ~ ^{r<N-i}^{l<j<i} ~ ^{N+l-i<r}^{l<j<N-r}- 



(3.22) 



10 



For example, for j = 0, we have: 

Vr = 2,...,Ar_i, Arii,0) = 5{r<^}. (3.23) 

Let r € {2, ..., A^ — 1} and j G {1, ..., N}. Let us prove that: 

A(^>j) = ~^{0<i~j<r-l}^{j<N-r} = ~^{j<i}^{i-j<r-l}^{j<N-r}- (3.24) 

Let us see the four following cases: 

• if J > i, then j > i — r. Moreover, ifA^ + 1 — f < r then j > N — r. In this case, 
Ar{i,j) =0-0-0 = 0. 

• ifi— j>r — 1 then r < i, j < i — r, j < i and j < N — r. Then, Ar{i,j) = 1 — 1 — = 
or Ar{i,j) = 1—0 — 1 = depending on r <N — i or N+1 — i<r. Finally, in this 
case, Ar{i,j) = 0. 

• if j > N — r then j > i — r. Moreover, if r < N — i then j > i. In this case, 
Ar{i,j) =0-0-0 = 0. 

• ^i j ^ i, i — j ^ r — l and j < N — r, then i — r<j. In this case, Ar{i,j) = 0-1-0 = —1 
or Ar{i,j) = — — 1 = —1 depending onr<A^ — zorA^ + 1 — i<r. Finally, in this 
case, Ar{i,j) = — 1. 

Consequently, (j3.24|) holds for any r € {2, ...,A^ — 1} and any j G {1,...,A^}. Finally, from 
<^2^ and (13723]) . we have: 

Vr = 2, ...,N - 1, \/j = 0, ...,iV, Ar{iJ) = S{j=o}S{r<i} - '5{0<i-i<r-l}'5{l<i<7V-r-}- (3.25) 

• Case i = or i = N. As Hq = 0, for any r = 2, ..., A^ — 1 and any j = 0, ..., A^, we define 
^r(0,j) = 0. As Hjy = Hj^^i, for any r = 2,...,N — 1 and any j = 0, ...,A^, we define 

Ar{N,j)=Ar{N-l,j). 

Hence, for any r = 2,...,N — 1, the elements (^Ar{i,j)) are defined for i = 0,...,N and 
j = 0,...,N. Then, we denote by Ar the matrix (^r(*)i))o<i xjv ^ -^A^+i- Finally, from 
(fSTTl) . we have proved that Al ( - aiF ■ a{G) + H) = where H = Y.^'^ «r-^r [F ■ a''(G)) . 
Finally, denoting Ai = —Idjsi+i G Mn+i-, we conclude the proof. D 

Remark 12. Using the discrete Leibniz formula given by ()3.14p . another choice in order to 
give a discrete version of (|2.16p from Lemma would be to apply the discrete version of the 
method used in Section \2.3[ Nevertheless, we would encounter many numerical difficulties. 
Firstly, such a method would imply the use of the operator A^ but this operator approximates 
the operator {d/dty only for tk with k > p. For the p first terms, we obtain in general a 
numerical blow up. Secondly, the use of operator A^ implies the use of h^^ . Hence, for a 
large enough p, we exceed the machine precision. 
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3.4 Numerical test 



We consider the classical bi-diniensional example (d = 2) of the quadratic Lagrangian L(x, v, t) ■ 
(x^ + f ^)/2 with [a, b] = [0, 1] and a = 1/2. Then, L is A"-invariant under the action of the 
rotations given by: 



{s,Xl,X2) 



cos(s) — sin(s) 
sin(s) cos(s) 



Xi 
X2 



The computation of (EL?) gives the following graphics: 



(3.26) 



We choose N = 600, Qo = (1, 2) and Qn = (2, 1). Let F and G denote: 

F = |^(0,Q) = (-Q2,Q1) and G = |^(Q, A^Q,t) = (A^Q^, A^Q2)_ (327) 




-0.15 
-0.151 
-0.152 
-0.153 
-0.154 
-0.155 
-0.156 
-0.157 
-0.158 
-0.159 

-0.16 



- constant of motion 



(a) Solution Q = (Q\Q'^) of pLg| 



(b) Quantity E^"/ Or A. {F ■ a'-(G)) 



Figure 1: Computation of (EL^) 



Figure [Ta| represents the discrete solutions of ( EL^ ) 



• and Figure [Tbl represents the explicit computable quantity Ylr=i OirAr[F ■ a'''{G)). As 
expected from Theorem [TOl this quantity is a constant of motion of ( jEL^p . 
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